Theoretical Investigation of Optical Conductivity in Ba(Fei_ x Co :r )2As 
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We report on theoretical calculations of the optical conductivity of Ba^ei-^Co^As, as obtained 
from density functional theory within the full potential LAPW method. A thorough comparison 
with experiment shows that we are able to reproduce most of the observed experimental features, in 
particular a magnetic peak located at about 0.2 eV which we ascribe to antiferromagnetic ordered 
magnetic stripes. We also predict a large in-plane anisotropy of this feature, which agrees very well 
with measurements on detwinned crystals. The effect of Co doping as well as the dependence of 
plasma frequency on the magnetic order is also investigated. 
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PACS numbers: 

The recent discovery of superconductivity in iron- 
based pnictide compounds 1 - has triggered intense experi- 
mental and theoretical research activity. A great deal of 
attention has been devoted to the relationship between 
superconductivity and magnetism. Similar to cuprates, 
the parent compounds exhibit antiferromagnetic spin- 
density-wave (SDW) order that disappears upon doping 
or pressure, giving rise to superconductivity. A possi- 
ble coexistence of magnetism and superconductivity has 
also been observed^, although this issue is still under 
debated. 

Despite the many similarities to cuprates, important 
differences characterize the pnictides; most relevant be- 
ing that their parent magnetic phase is not Mott insu- 
lating but rather metallic in nature. The preponderance 
of magnetic properties of FeAs compounds can be de- 
scribed by local/semilocal approximations to exchange 
correlation functionals within density functional theory 
(DFT), however the debate on the nature of magnetism 
(itinerant or localized) in this class of materials is still an 
open question^. The degree and type of correlation in 
these materials is also under debate, due mostly to the 
ARPES data which suggests a strong renormalization of 
the Kohn-Sham energy bands^. 

The pairing mechanism in pnictides is not yet uni- 
versally agreed upon - but the possibility of standard 
phononic superconductivity has now been ruled ou t 15 ' 16 . 
It has been suggested that the superconductivity in these 
materials could be mediated by magnetic excitations cou- 
pling electron and hole pockets of the Fermi surface and 
favouring s-wave order parameters, with opposite sign on 
different sheets of the Fermi surface (s± coupling) . 

What is clear is that the main physics in FeAs materi- 
als is controlled by a subtle interplay of magnetism and 
Fermi surface topology. One of the tools that can be used 
to probe these properties is optical measurement. In this 
work we report on our results of the optical properties 
of pure and Co-doped BaFe2As2 using DFT. A detailed 
comparison of this work with infrared optical conductiv- 



ity experiments is also performed, and provides interpre- 
tation of most of the main features of the experimentally 
observed spectra^ - —. 



I. METHODOLOGY 

Our calculations were performed using the highly ac- 
curate all-electron full potential linearized augmented 
planewave method (LAPW) 23 as implemented within the 
Elk 24 code. 

All calculations presented here have been performed 
using the experimental crystal structure for undoped 
BaFe2As2^, and not the low temperature orthorhom- 
bic distortion. The main spectral features for the doped 
BaFe2As2 have been checked to be independent of this 
choice. 

The optical conductivity tensor ct(w) is related to the 
dielectric function b}* 2 ^ (atomic units are used through- 
out): 



47T2 

Cc^M = $aj3 H ^Mi 



(1) 



for Cartesian components a, ft. In practice, a is calcu- 
lated according to^l 



O 2-^ cv 



VCK 



n« c (k)n^(k) , (n« c (k)n^(k))* 



IT] 



u - e vc + 



4tt 



uj - 



irj 



ri] 
(2) 



where the first and second terms are the interband and 
the intraband contributions respectively; e vc = e v —e c are 
the differences between valence and conduction eigenval- 
ues; ^ are /c-point weights over the Brillouin zone; f2 is 
the unit cell volume; and 
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are the momentum matrix elements. The imaginary term 
ir\ mimics a broadening of the spectrum owing to a finite 
lifetime, and a value of 0.1 eV is used for this. The plasma 
frequency, which is also a tensor quantity, is given by 
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with ep being the Fermi energy. We find the dielectric 
function and plasma frequency to be sufficiently well con- 
verged with a mesh of 16 x 16 x 4 fc-points. 



II. RESULTS AND DISCUSSION 



A. Undoped compound 



1, Ground state properties 



We investigated the non-magnetic (NM) phase and 
various possible magnetic orderings. The most relevant is 
the striped phase (Str.), where on each Fe plane stripes 
of up spins are alternated with stripes of down spins. 
We also studied the ferromagnetic (FM), checkerboard 
(CB) and non-collinear antiferromagnetic (2K) order- 
ings. In the latter case, spins are rotated 90 degrees 
between adjacent Fe sites. The corresponding total en- 
ergies, relative to the Str. phase, have been collected in 
Table HI In agreement with experiments and other DFT 
calculations^, we find the striped phase to be the most 
stable. 





NM 


CB 


Str. 


FM 


2K 


E (meV) 


120.0 


70.1 





119.5 


33.1 


M(Fe) (/i S ) 





1.428 


1.490 


0.31 


1.33 


LU P x, y 


2.55 


1.62 


0.702/0.91?/ 


1.958t/l-9634. 


1.44 


lo p z (eV) 


0.79 


1.54 


1.03 


0.60170.874, 


2.349 


vf x,y 


13.56 


8.94 


7.77z/10.07i/ 12.611723.054. 


9.168 


v F z (10 4 m/s) 


4.22 


8.52 


11.45 


5.961710.264- 


14.915 



TABLE I: Undoped BaFe2As2 data for different magnetic or- 
derings. E is the total energy per formula unit relative to the 
Str. phase; M(Fe) is the magnetic moment integrated in the 
Fe muffin tin; top and vj are the plasma frequency and Fermi 
velocity respectively. In-plane components are labeled by x, y 
and the out of plane component by z. For ferromagnetic or- 
dering the symbol "\ indicates the majority spin channel. In 
the Str. phase the spin stripes are oriented along the x axis. 



As in the case of LaOFeAsii, the 2K structure is close 
in energy to that of the striped, while the FM solution 
is almost degenerate with the non-magnetic one. The 
magnetic moment, integrated in the Fe muffin tin, is 
larger than 1 /jb in all three antiferromagnetic struc- 
tures, and is equal to 1.49 /i^ in the ground state Str. 
phase (Tab. These theoretical moments are much 
larger than the experimental values measured by means 



of neutron diffraction (« 0.9/is— iS 9 .— ) and Mossbauer 
absorption(0.4 /j,b^)- 



2. Plasma frequencies 

We first focus on the plasma frequency tensor (w p ), 
that determines the low frequency part of the spectrum. 
Owing to lattice symmetry, oj p is diagonal with respect to 
the tetragonal axis, and we will indicate the xx, yy and zz 
components simply as x, y and z. In all the investigated 
magnetic structures u> p has only two independent com- 
ponents, apart from the Str. phase which has three. For 
the FM case we define separated plasma frequencies for 
majority and minority spin channels. The corresponding 
results are listed in Table fl] 

The non magnetic phase yield a value to p of 2.55 eV. 
Magnetic order strongly affects the value of the plasma 
frequency in that magnetism opens a partial gap at Ep 
thereby reducing both the Fermi velocity and the density 
of states at the Fermi energy 3 - 3 - and resulting in a decrease 
of uj p as compared to the non magnetic case. The antifer- 
romagnetic phases have the most isotropic u> PtZ /uj p ^ x ra- 
tio, while the NM structure is the most anisotropic. Our 
non-magnetic plasma frequencies are in agreement with 
previous calculations 3 -! - — . Magnetic plasma frequencies 
in the Str. phase have been reported before in two recent 
works 3 §^£ and we find a reasonable agreement with our 
results. Slightly smaller and more isotropic values are 
reported in the past perhaps because of the larger value 
of the Fe moment (1.64 /is and 1.98 /xb). 

The experimental in-plane plasma frequency, extracted 
from optical conductivity measurements, has been esti- 
mated in Ref. HI to be w p = 1.6 eV at room tempera- 
ture and 0.58 at 10 K. The low temperature limit is in 
good agreement with our estimate for the striped phase 
(our zero temperature ground state). On the other hand, 
the room temperature value, measured above the Neel 
temperature (TV), should be compared to our estimate 
in the non-magnetic phase (2.55 eV); the experimental 
value turns out to be much smaller. A temperature ef- 
fect, not included in our calculations, could partially jus- 
tify this disagreement. However a strong dependence of 
the plasma frequency on temperature is not expected: in 
case of a similar system, BaNi2As2, the plasma frequency 
drops by only 5% of its value between 300 K and 150 K 
(above the structural/magnetic phase transition) 39 . 

For the non-magnetic phase, a disagreement between 
LDA calculated plasma frequencies and experimental es- 
timation from optical conductivity has been already re- 
ported, and been interpreted as an indication of correla- 
tion effect o 34 i 35 ' 40 i 41 . However, assuming that the main 
reason for the discrepancy is correlation, our results in- 
dicate that renormalization effects are much stronger in 
the normal state than in the magnetic one. If we con- 
sider the ratio between the theoretical and experimental 
w^, asa measure of the renormalization effect induced by 
correlations, as was done in Ref. 1341 . then quite unexpect- 
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FIG. 1: (color online) Optical properties for non-magnetic 
BaFe2As2. In the top panel Real (black) and Imaginary (red) 
components of the dielectric function. Dashed and full lines 
refer, respectively to intraband and interband contributions as 
in eq. [2] In the two central panels, the optical conductivity 
tj x , y refers to the in plane and a z to the out-of-plane compo- 
nent. Dashed lines are the intraband contributions. The joint 
density of states is resolved into orbital contributions in the 
lowest panel. 



edly we find the non-magnetic phase of BaFe2As2 to be 
a strongly correlated material, while magnetic BaFe2As2 
appears to be a conventional one. 

There has been disagreement in the literature on 
whether this implies that the dynamical spin correlation 
effects present in the metallic phase are frozen out in the 
magnetic one, or that the effect can not be ascribed to 
correlation effects at all, but is simply due to the inad- 
equacy of the LSDA exchange-correlation functional (as 
suggested by Mazin in Ref. 41). The veracity of either 
of these scenarios cannot be determined by the present 
calculations 



3. Optical Conductivity 

Calculated optical conductivity for the non-magnetic 
phase is shown in the middle panels of Fig. [TJ Strong 
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FIG. 2: (color online) In plane optical conductivity for various 
magnetic configurations of BaFe2As2. Full and dashed lines 
give respectively real and imaginary parts of a. All the data 
refer to the interband components only, apart from the FM 
case where the Drude contribution is also reported as a dotted 
line. 



anisotropy between the in-plane (a xy ) and out of plane 
(ct z ) conductivity is immediately apparent. In particular, 
the out-of-plane conductivity has a more pronounced fre- 
quency dependence than its in-plane counterpart. This 
anisotropy underscored by the appearance of a plasmonic 
peak, present only in a z . At this point it is worth noting 
that the value of the plasma frequency in the out-of-plane 
case is much smaller than in conventional metals, the rea- 
son for this lies in the low value of the Fermi velocity in 
the z direction. 

The origin of optical conductivity in the NM phase can 
be determined by looking at the joint density of states 
(JDOS), which is the real part of the optical conductivity 
with equal weight for all optical transitions (momentum 
matrix elements set to 1). This JDOS can be further 
resolved into various orbital contributions as follows: 

JDOS(£, l)= J2 6 ( £ " k _ £ »' k _ £ ) ™kP„'k(0 (5) 

where n are band indices and P n k(l) is the projection 
of the Kohn-Sham state on the orbital channel I. This 
orbital projected JDOS is shown in the lowest panel of 
Fig. [TJ Clearly, all the low energy features originate from 
the Fe (i-states that dominate near the Fermi Energy. 

The optical conductivity in various magnetic phases is 
shown in Fig. [5] The main effect of magnetism is an 
increase of conductivity in the 0.1 - 1.0 eV range. The 
reason for this is that the reduced symmetry makes more 
states available (folding of the Brillouin zone due to the 
antiferromagnetic symmetry breaking), and new inter- 
band transitions are allowed in the low energy range. As 
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A rX M rZR AM striped phase of BaFe2As2. Defined as the difference between 

conductivity in the off-stripe direction and stripe direction. 
Different curves correspond to a different magnetic moment at 

FIG. 3: (color online) Band structure near the Fermi energy the Fe site (see legend for the moment value in atomic units). 

(tetragonal unit cell). Upper panel: comparison between non- Drude contribution is included. 



magnetic and magnetic stripe ordered phase. Lower panel: in 
the striped phase, comparison between the undoped (x — 0.0) 
and Co-doped (x = 0.03) magnetic structure. 



we will see in the following, the new structures merge in 
the Drude intraband contribution, as soon as the mag- 
netic order disappears. This is anticipated by the be- 
havior of a in the Ferromagnetic phase, where the or- 
dered moment is low. In this phase a broad shoulder 
at about 0.1 eV is almost merged with the Drude peak 
(two dotted lines corresponding to the f and \. channels). 
This enhancement of cr(w) in the 0.1 - 0.3 eV region, in- 
duced by AF ordering, is an important feature of our 
results. This has been observed experimentally by differ- 
ent groups. We identify the peak at 1000 cm -1 (0.124 
meV) reported in Refs. H7lll8ll20ll2ll as a result of the 
AF ordering which introduces additional low frequency 
transitions. It is interesting to notice that in the striped 
phase it appears only for field polarization perpendicular 
to the stripes (let's call it off-stripe). No similar effect 
exists when the electric field is in the direction of the 
stripes where the magnetic ordering is ferromagnetic. 

In order to understand this we plot in Fig. [3] the en- 
ergy bands dispersion along the main lines of the Bril- 
louin zone, including the two orthogonal in-plane direc- 
tions both in the top and basal planes. Compared to 
the usual representation of bands in this compound, we 
have here a folding of the Brillouin zone, with the M 
point electron pockets folded back at T. We highlight in 
particular a Fe d— band located at about —0.1 eV at T, 
strongly dispersed in the direction of the stripe (r — Y) 
and flat in the off-stripe direction. The presence of this 
band offers a much larger phase-space availability for ver- 
tical transitions, compared to the non-magnetic case and 
to the orthogonal direction. If this interpretation is cor- 
rect, then the experimental peak can be fully understood 
in terms of single particle properties. 

At higher energies, around 8 eV, a peak mostly due to 
Ba /-states appears in the imaginary part of er, inducing 



a corresponding structure at about 6 eV in its real part. 
These high energy features are not affected by the mag- 
netic ordering, which induces important major changes 
only up to about 3 eV. 

The effect of magnetic ordering can be singled out by 
defining the conductivity anisotropy Aer(oj) = o-y(us) — 
<j x (uj), which is plotted in Fig. H) This anisotropy is 
strongly dependent on the magnetic moment, which in 
turn is strongly overestimated by standard DFT-LDA. 
In order to partially account for this effect, we perform 
fixed-spin-moment DFT calculations and report Act as a 
function of the Fe magnetic moment varying from 1.5 [Ib 
(the value obtained in our ground state calculations) to 
zero. 

We see that for the calculated ground state, the low 
frequency part of the spectrum has an anisotropy peak 
centered around 0.25 eV (« 2000 cm -1 ). The spectral 
weight of this peak is compensated by depletion regions 
at high energy (see for example the negative peak around 
1.5 eV). For smaller values of the magnetic moment this 
peak grows in intensity and also shifts toward lower fre- 
quencies. According to the above interpretation of the 
peak of tr, this can be understood in terms of the cre- 
ation of a quasi-gap (depletion in the density of states) 
around the Fermi level Ep : a smaller moment implies a 
smaller quasi-gap, with a more peaked shape anisotropic 
structure. The same happens to the high energy negative 
structures. 

Our theoretical spectrum and its anisotropy can be 
compared with the experimental results^. The presence 
of anisotropy peak, magnetic in origin, is confirmed by 
experiment, both in its energy position and in its general 
shape, with a depletion of states around and above 500 
meV (4000 cm -1 ). To have a consistent comparison, we 
should consider our values calculated at the experimen- 
tal magnetic moment. This is not easy as there is no 
universal agreement between different experiments: neu- 
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FIG. 5: (color online) Top panel, Experimental optical con- 
ductivity of detwinned Ba(Fei_ 2 ,Co a: )2As measured with the 
field along the two in plane Cartesian directions (red and blue 
dashed lines, compared with theoretical prediction in the mag- 
netic Str. phase (full lines). Bottom panel, experimental op- 
tical conductivity anisotropy (a y — a x ), at 0.025 of Co doping 
and 10 K, compared with the theoretical result obtained for 
a reduced magnetic moment of 0.75 /is, close to experimental 
estimations (see text for details). 



tron diffraction and Mossbauer measurements give mo- 
ments around 1 and 0.4 fis respectivel y 19 ' 29 " —. In light 
of this, our comparison can only be semi-quantitative. 
Still, using the results for fip e = 0.75 — 1/ab, the peak 
at 1000-2000 cm' 1 (120 - 250 meV) and the depletion 
of states after 3000-4000 cm" 1 (350 - 500 meV) found in 
experiments are in fact well reproduced by theory. We 
also notice that in experiments the peak position does 
not change with temperature, while in our case it does 
change with the magnitude of the fixed moment. While 
the latter finding is clearly a consequence of the evolution 
of the band structure, the experimental behavior seems 
to point to a disorder of atomic moments, rather than to 
their reduced magnitude. 

Bearing this in mind, we can make a detailed com- 
parison between our results and experiments based on 
the anisotropy in the optical conductivity recently ob- 
served on detwinned crystals^ 8 -. In Fig.[3]we plot the ex- 
perimental anisotropy together with theoretical results 
(Drude contributions are not included in both curves). 



We clearly see a fairly good agreement in the general 
shape. The main interband peak starts with a slightly 
broad shoulder at lower energy in experiments than the 
theoretical calculations. Most importantly, by compar- 
ing the xx and yy components of a we see very simi- 
lar features in the magnetically induced anisotropy. The 
peak is found at higher energy in our theoretical results; 
however, as seen in the previous discussion (see also Fig. 
3]), this peak moves to lower energies as the Fe magnetic 
moment reduces. If we force the magnetic moment to be 
close to the experimental value (see lower panel in Fig. [5]) 
the comparison improves significantly, with a very good 
agreement between theory and experiment. 



B. Co doping 

A fuller understanding of the optical properties of this 
system can be achieved by studying the doped system. 
It is well established experimentally that electron or hole 
doping of the BaFe2As2, rapidly destabilizes the mag- 
netic phase. In the case of Ba(Fei_ a: Co a: )2As at 6% Co 
concentration magnetism disappears and the compound 
becomes a superconductor. A possible coexistence be- 
tween the two phases has been report e d 5 ' 8 ' 42 , possibly 
justified by nano-scale inhomogeneity in the dopants. In 
this work we assume that the magnetic phase is destroyed 
at 6% of homogeneous Co doping, and that below this 
doping concentration superconductivity is absent. 

Keeping the crystal structure fixed, and using the vir- 
tual crystal approximation (VCA), this phase transition 
is not properly reproduced by DFT LSDA calculations 33 . 
The Fe moment seems to be only weakly reduced by elec- 
tron doping (by about 20% for x=0.06). The strong sensi- 
tivity of magnetism to the As positions and the unavoid- 
able impossibility of describing disorder effects could be 
all valid reasons for this disagreement with experiments. 
This disagreement could again be due to correlation ef- 
fects beyond LSDA. We can simulate the experimental 
situation (and in particular the Co induced destabiliza- 
tion of magnetism) by doing fixed spin moment calcula- 
tions, as before. At each doping (achieved using VCA) 
we fix the magnetic moment by scaling it according to the 
experimental magnetic critical temperature Tn*&. This 
constraint is likely to correct the deficiency of the con- 
ventional LSDA calculation and will mimic the two main 
Co induced effects: a localized electron doping on the Fe 
layer, and a magnetic to non-magnetic phase transition 
induced by this. 

The calculated optical conductivity in the Str. phase is 
plotted in Fig. [6] When doping reduces the Fe magnetic 
moment the system becomes more metallic - an increase 
of the plasma frequency is accompanied by a large trans- 
fer of spectral weight from the structures around 0.8 eV 
to the Drude peak. The region of suppressed conductivity 
around 70 meV, is filled as the magnetic gap closes. The 
observed trend is similar to the one reported in experi- 
ments (see for example the work of Nakajima et al. Hi] 



6 




0.1 1 

co(eV) 



FIG. 6: (color online) In plane optical conductivity 
Ba(Fei_ x Co :E )2As, in the Str. phase. The x direction cor- 
responding to the direction of the stripes. The x and y direc- 
tions are averaged in the bottom panel, as they should come 
in an experimental measurement for twinned samples. Drude 
contribution is included. 



and Lucarelli et al. fl7t ). The main difference between our 
calculations and the experimental work (see for example 
Fig. [4] on Ref. .21), is in the evolution of the spectra: we 
observe a transfer of spectral weight from one peak in the 
high energy part of the spectrum to the Drude compo- 
nent, while in experiments the spectral weight transfers 
smoothly, without spectral peaks moving sharply as a 
function of x. 

We believe that this difference can be also understood 
in terms of the intrinsic inhomogeneity of the doping in 
samples, present even in a macroscopically homogeneous 
Co distribution. Experimental indications of this kind 
of inhomogeneity have been reported and discussed in 
several papers^. Therefore, experimentally each spec- 
trum would correspond, to some extent, to an average 



of a range of theoretical dopings, that are intrinsically 
homogeneous in the virtual crystal approximation. 



III. CONCLUSIONS 

In this work we report the calculated optical conductiv- 
ity for BaFe2As2, one of the most studied members of the 
iron arsenide family in order to understand which exper- 
imental features can be reproduced by theoretical meth- 
ods based on DFT. We confirm the fact that the plasma 
frequency in the non-magnetic phase is much higher than 
the experimentally observed values, indicating a poor de- 
scription of the non-magnetic state within LSDA. On the 
other hand, in the magnetic striped phase the experimen- 
tal estimate of the in plane plasma frequency (0.58 eV) is 
not far from our values (0.7 and 0.91 eV). If one ascribes 
the disagreement between theoretical and experimental 
results to correlation effects, as proposed by Qazilbash 
et al, this would mean that the virtual spin fluctuations 
freeze in at the onset of static magnetic ordering. 

Studying the frequency dependence of the optical con- 
ductivity, we observe that a low energy peak appears 
whenever the polarization of the incident light is directed 
along an antiferromagnetically oriented spin configura- 
tion, as in the CB or Str. phase along the y direction. 
This magnetic peak has been observed experimentally by 
various groups. In this work we are able to explain the 
experimentally observed large anisotropy of this peaki^. 

We also reported the effect of Co doping, by combining 
fixed moment calculations (to reproduce the experimen- 
tal Fe moment) and the virtual crystal approximation 
and were able to reproduce the spectral weight transfer 
of the optical conductivity from the high to the low en- 
ergy parts of the spectrum, as a function of such doping. 
This effect can be ascribed to the reconstruction of the 
electronic structure around the Fermi energy because of 
the magnetic ordering. 
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